Different ocular dominance map formation influenced by 
orientation preference columns in visual cortices 
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In animal experiments, the observed orientation preference (OP) and ocular dominance (OD) 
columns in the visual cortex of the brain show various pattern types. Here, we show that the 
different visual map formations in various species are due to the crossover behavior in anisotropic 
systems composed of orientational and scalar components such as easy-plane Heisenberg models. 
We predict the transition boundary between different pattern types with the anisotropy as a main 
bifurcation parameter, which is consistent with experimental observations. 
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The highly ordered structure in the mammalian vi- 
sual cortex has attracted much attention from theoreti- 
cal neurobiologists and has been thoroughly studied with 
the expectation of providing the basis for neural dynam- 
ics and computational models. Though most models of 
the visual map formation are based on common postu- 
lates, such as Hebbian synapses, connections or competi- 
tions between neighboring neurons and synaptic normal- 
ization, there are quite a number of successful models 
with unique mechanisms 0, . The Hamiltonian models 
with spin variables were proposed for the visual map for- 
mation with a striking analogy with the physical systems, 
such as magnetism 0- Recently, the characteristics 
of visual maps are systematically explored through the 
statistical properties of vortices in magnetism ||. The 
spin-like Hamiltonian models represent essential ingredi- 
ents of neural interactions in the visual map formation 
without paying much attention to particular neural con- 
trol mechanisms and can be shown to exhibit common 
statistical properties of the vortex formation as in other 
development models [(J. 

As vast experimental data on visual maps are accu- 
mulated, the various patterns in visual cortices of dif- 
ferent animals have drawn much interest with the ex- 
pectation of testing different neural models experimen- 
tally 0, H H E O 13 El The observed visual pat- 
terns can be classified as at least three different types. In 
macaque monkeys, the OD columns form parallel bands 
of regular spacing with relatively few branching points 
that are mainly oriented perpendicular to area bound- 
aries Q. The degree of OD segregation is strong and 
the typical average spacing of OD, Aod, is larger than 
that of the OP columns, Aop (Aod > Aop) M- ^ n cats 
and ferrets, OD columns form an array of beaded bands 
exhibiting only a weak tendency of elongation perpen- 
dicular to area boundaries 0, 0, 0] . The degree of OD 
segregation is intermediate with Aod < Aop [l^. In 
the case of tree shrews, the OD segregation is very weak 
or absent, while more stripe- like patterns are observed 
in extensive OP columnar regions with low densities of 



orientation centers |13| . 

Some experimental works show that the OP and OD 
patterns are structurally correlated. The singular points, 
so-called pinwheels, in OP columns tend to align with 
the centers of OD bands and the iso-orientation contours 
intersect the borders of OD bands at right or steep an- 
gles |8(. Such correlations between two columns occur 
due to the normalization of the synapse strength [l4| . 
The average response of the orientation selectivity van- 
ishes at singularities, so that the strong OD components 
can develop near pinwheels. The orthogonality between 
contours of the columnar patterns can be explained intu- 
itively through the isotropic Heisenberg model 5]. The 
different visual pattern types originate from both the 
synaptic normalization and the anisotropic interaction 
between two columns. The parameter A, a measure of 
the anisotropy between the columns, turns out to be cru- 
cial in determining not only the pinwheel stability but 
also OD pattern types and OP pattern regularity. We 
estimate the anisotropy for each species from the typi- 
cal spacing of OP and OD patterns from experiments. 
We calculate the threshold of anisotropy between differ- 
ent OD pattern types and construct the phase diagram, 
which are consistent with experimental observations. 

One of reasons why so many different neural models 
lead to the successful formation of visual maps is that the 
typical characteristics of self-organizing maps are deter- 
mined by the topology of lattice and feature space rather 
than the detailed cortical modification rules. The sta- 
tistical properties of emergent cortical maps can be pre- 
dicted using only minimal mathematical constraints such 
as the symmetry. In the fiber bundle map (FBM) repre- 
sentation method, the phases of feature components are 
described by the continuous group corresponding to the 
manifold in the feature space Because OP columns 

have 0(2) (or U(l)) symmetry, the energy functions for 
the OP map formation should be invariant under gauge 
transformations and take a general form 



dr |-|(V - iA)^ OP \ 2 - — \^op\ 2 + • • • J 
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for the oricntational feature components -0op(r) = 
(q(r)cos20(r),g(r)sin2< ? (>(r)) (or q(r)e 2 '^ (r) ) with the 
preferred angle 4>(r) and the degree of preference q(r) at 
the cortical location r. This can be obtained also from 
other map formation models after a continuum approx- 
imation and can explain most typical properties of OP 
patterns found in experiments and simulations 

Biologically the cortical modules are composed of sev- 
eral layers and the different feature components are found 
from the selective response properties of corresponding 
layers. Each layer can have a different synaptic strength, 
so that the energy for the combined OP and OD map 
formation takes an anisotropic form 

E = |dr|^|(V-zA p)^ p| 2 -^|^ p| 2 (2) 

+ ^|(V - zA OZ5 )^ OD | 2 - ^piT/WI 2 } 

including the scalar component representation of OD 
ipOD- The restriction imposed by the normalization for 
the total feature vector i/>(r) = (?/>op( r ), ^oz)( r )) is that 

\ij(r)\ 2 = \^J p(r)\ 2 + I^odWI 2 = const (3) 

for all r. The normalization and the equilibrium (V 2 i/' AI ~ 
for fjb = OP and OD) conditions lead to the or- 
thogonality between OP and OD contour lines, that 
is \7ijJop ■ V^od ~ 0. Now, we define a parameter 
A = mQ D /rriQp describing the anisotropy between two 
columns. Then Eq.|0 has two different solutions de- 
pending on the anisotropy A. If A < 1 (or A > 1), it 
has a ground state at the energy density — TOq P /2 (or 
—m 2 OD /2) with weak OD segregations (tpo D ) = (or 
strong OD segregations (V'od) = l^l 2 )- The previous 
studies in magnetism have shown that there exists an- 
other threshold A c when features are composed of ori- 
entational and scalar components. The statistical prop- 
erties of visual maps are closely related to those of the 
classical two-dimensional anisotropic Heisenberg model 
(CTDAHM), described by the Hamiltonian 

H = -J^SfSf + SfS* + XSfSj), (4) 

for J > 0. The possibility of two vortex types, named 
in-plane and out-of-plane, in the CTDAHM was first dis- 
cussed by Takeno and Homma [15J. Different values of 
the critical anisotropy A c , above which out-of-plane vor- 
tices become stable, were found to exist and depend on 
the lattice type (A c w 0.72, 0.82, 0.62 for square, honey- 
comb, and triangular lattices, respectively) flif . 

In this Letter, the anisotropic Heisenberg model with 
distance-dependent interactions is used for simulations : 

H = - { D op(rij)(SfS* + S?S») + Do D {nj)S?S*}(5) 
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(c) k = 0.3, A = 0.63 (d) k = 0.3, A = 0.67 



FIG. 1: Development of the scalar peak near the orien- 
tational singularity. The preferred angle <j> is denoted in 8 
colors and the ocular dominance S z is denoted in black and 
white wired frames, where the black denotes the region with 
Sz > and the white, otherwise. Simulations are done for 
the Hamiltonian in Eq.lJKJ with oq D — oq P = 6, koD = kop 
(A = eod/eop)- (60 x 60 lattice and a non-periodic boundary 
condition) 

for the preferred angle & = (1/2) tan" 1 ^ /Sf ), the OD 
Sf with normalization to a unit modulus (|Si| 2 = 1). 
This satisfies also the energy form in Eq.|J2J in a con- 
tinuum approximation. Other development models, such 
as the elastic-net model or Kohonen's SOFM algorithm , 
can be rewritten into the spin- like Hamiltonian model ((J . 
Usually _D(x,y) = -D(|x — y|) is positive in short-range 
and negative in long-range, that is, of the Mexican hat 
type. The typical parameters in Eq.Q or and the 
exact anisotropy A are determined by the actual form 
of D(r). In general, the interaction strength within each 
column is proportional to the activity strength e M , so that 
A oc Sod/^oPi and to the number of interacting pairs, so 
that A oc Cq^/cqp in two-dimension for the cooperation 
range a^. We take the form of the interaction function 

D»{r) =e»(l- kp^j exp (-r 2 /2 ( x 2 ) (6) 

for an inhibitory strength Single vortex simulations 
in Fig. show the structure of two vortex types. When 
A is above the threshold A c , a peak in the OD map starts 
to develop near a singular point in the OP map (Fig. 
and^i) . When there exists a weak inhibitory activity, we 
have ferromagnetic solutions with both a scalar peak and 
and an orientational singularity (Fig. \3>), which are the 
typical features of in-plane and out-of-plane vortices ob- 
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FIG. 2: (a) Different OD pattern formations depending on 
the anisotropy A. The gray color in boxes denotes the weak 
segregation region (\S Z \ < 0.01). Parameters are u p — 
a OD = 6, kop — koD — 1.0 (A = eod/eop)- Maps are 
generated using the Metropolis algorithm at zero tempera- 
ture, with a non-periodic boundary condition and an initially 
random state in the 70 x 70 lattice. 



served in the CTDAHM. Similar results can be obtained 
in different map formations models, if they are modified 
to have an anisotropy between columns. 

Fig. [21 shows the emergence of three different pattern 
types in OD columns; Type (A) for A > 1, parallel bands 
with strong OD segregation, Type (B) for A c < A < 
1, beaded bands with intermediate OD segregation, and 
Type (C) for A < A c , the absence of OD segregation. 
Types (A), (B), and (C) in Fig.|21correspond to observed 
OD patterns in macaque monkeys, cats, and tree shrews, 
respectively. We measure the degree of OD segregation 
in terms of the standard deviation of S z at a given time 
that is averaged over 10 trial evolutions for a given A. 
We find that the OD segregation strength increases as 
A increases. The in-plane vortices are unstable under 
small thermal fluctuations for A < 1. The anisotropy 
A is related to the pinwheel stability and connected to 
the measured pinwheel density for different animals 5, 
Il7j . Beaded bands are not permanent structure as well 
and disappear as pinwhccls annihilate during the map 
development process. 

For A > 1, which corresponds to the Ising model re- 
gion, the OD columns play a dominant role in the over- 
all map formations. Parallel bands appear in a whole 
range of the OD map and make a right or steep an- 
gle with area boundaries, which are the typical prop- 
erties of OD maps generated without OP columns. In 
the OP columns, singular points and periodic patterns 
still exist, but are much less regular. The distribution 
of iso-orientation contours is not uniform with locally 
dense or sparse regions (see Fig. The rotational 

plane defined from the gradient directions in the feature 
space helps to understand this point. For A = 1, the 
map formation process is balanced between two station- 
ary solutions with the normal vector of the rotational 



plane slanted at \n x \ = \h y \ = \h z \ = l/\/3, so that 
AS Z = V 1 - \n z \ 2 /V2 = 1/V3 0.577, as seen in 
Fig. |21 For A > 1 (or A < 1), the rotational plane lies 
more vertically (or horizontally), so that AS Z approaches 
1 (or 0) as the map formation progresses. If the rotational 
plane becomes too slanted, the circular trajectories are 
projected to elongated ellipses in the S x — S y plane, lead- 
ing to irregularities in <f> coordinates. 

In the real brain, it is not certain which factor among 
various parameters such as the activity rate e M , the co- 
operation range ct m and the inhibitory strength k^, is a 
major contributor to the anisotropy between OP and OD 
columns. We can obtain some insights on this by com- 
paring the typical spacings Aop and Aod measured in 
animal experiments. Various models predict that the typ- 
ical spacing A M increases in proportion to the cooperation 
range (A„ oc Ou) and decreases for stronger inhibitory 

strength £^H[ll[l^]. If 0P and 0D columns are de- 
ferent only in the activity rate e^, the typical spacings 
are the same, that is Aod — Aqp for all A. If they differ 
only in fc M , the typical spacing of OD bands is usually 
smaller than that of OP slabs when the strong OD seg- 
regation occurs (Aod < Aop for A > 1). However, if 
they differ only in er M , the typical spacing of OD bands 
is larger than that of OP slabs when strong OD segre- 
gation occurs (Aod > Aop for A > 1). Among these 
possibilities, we find that the last case is consistent with 
the animal data in cats JJ, [HEl EH13 , ferrets [THI^ 
and macaque monkeys [8( . This leads us to use the coop- 
eration ranges between columns as the major bifurcation 
parameter for investigating different visual pattern types 
in animals. Fig. shows the results of simulations using 
the typical spacings Aop and Aod measured in animal 
experiments. For Aod < Aop, more beaded bands in 
OD columns emerge near the pinwheel centers as in the 
visual maps in cats (Fig. 0i), whereas for Aod > Aop, 
parallel bands in OD columns emerge in a whole region 
as in macaque monkeys (Fig.©). 




(a) (b) 

FIG. 3: Simulations of OD and OP maps for (a) Aod < Aod 
(a OD = 6.7, a 2 OP = 10.0) and (b) A OD > A OP (a OD = 13.8, 
a OP = 10.0). Parameters are eod = sop and kop = koD = 
0.7, so A = Ood/^op = A OD /A OP (100 x 100 lattice with a 
non-periodic boundary condition). 



4 



Parallel Bands 



Macaque monkey 



Beaded 
Bands. 




(Tree Shrew) 



Segregation Absent 



FIG. 4: The phase diagram of OD patterns according to the 
cooperation ranges, aoD and o~op, where dotted lines denote 
the experimental data for various animals. 



The exact calculation of the critical anisotropy A c , 
above which the out-of-plane vortices appear, is also an 
important problem in magnetism |2^. Unfortunately, 
any continuum theory fails in calculation of A c because 
of the singularity near the vortex core [25j. The de- 
velopment of out-of-plane components depends on not 
only the exact form of the interaction function D(r) but 
also other factors such as lattice types, vortex distribu- 
tions, etc. The boundary between the "OD segregation 
absent" and "Beaded bands" regions cannot be distin- 
guished sharply in both experiments and simulations be- 
cause the development of OD components is weak near 
this boundary. The analytic and simulational methods in 
the CTDAHM can be extended to calculate the critical 
anisotropy precisely for various cases j2(|. We find that 
the stability boundary of the out-of-plane vortex takes 
the form (tod — otaop + (3. Since the computed value of 
(3 is small, the region for the beaded band pattern can 
be described by a < <jod/vop < 1, where a ~ 0.6 for 
A = &od/ ij op- This agrees well with the experimental 
data, where Aod/^op ~ 0.82 for cats and ferrets. Fig.^ 
shows the phase diagram of different OD pattern forma- 
tions in the parameter of cooperation ranges, a op an d 
(Jod- The experimentally derived curves for macaque 
monkeys, cats, and ferrets lie well within regions of our 
predictions. 

In this Letter, we showed the influence of the correla- 
tion between OP and OD columns on the overall visual 
map formation, and explained the observed experimental 
data successfully. Moreover, we also showed that the typ- 
ical characteristics of emergent cortical maps can be de- 
scribed well by the solutions from the general field theory 



with minimal constraints. Our study of the visual map 
formation suggests another possibility in modeling phys- 
ical neuronal dynamics with physical models. In spite of 
the complex circuitry and nonlinear dynamics in neural 
systems, the neural behavior at high levels may obey sim- 
ple and general rules, which can be related to the noble 
theory of statistical physics. We expect to find more ap- 
plications of physics-based methodology in interpreting 
diverse neural phenomena and constructing noble com- 
putational architecture. 

This work was supported by the Ministry of Science 
and Technology and the Ministry of Education. 
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